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order in SCET. We present the factorized result for the decay amplitudes in rare radiative B decays 
at leading order in SCET, and at next-to- leading order in a s . 
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I. INTRODUCTION 



The formalism of soft-collinear effective theory (SCET) [1-6] can be applied to high- 
energy processes in which massless quarks are emitted with large energy. In this case a light 
quark with large energy has three different scales Q, Q\ and Q\ 2 , where Q is the large 
energy scale and A ~ P±/Q is a small parameter. Here p± is a typical transverse momentum 
scale. SCET has been applied to many high-energy processes in which final-state particles 
include energetic light particles. For example, the hard contribution to the heavy-to-light 
form factor [1, 2, 5, 7, 8], the pion form factor and high-energy scattering [9], and radative 
leptonic B decays B — > ^ev [10, 11] have been studied. Recently it has been applied to 
nonleptonic B decays into two light mesons by the authors [12, 13] and nonleptonic decays 
are shown to factorize at leading order in SCET and to all orders in a s . 

In this paper we apply SCET to rare radiative B decays B — > V7, where V can be a 
p meson or a K* meson. In this case, employing a two-step matching offers a convenient 
procedure to obtain the matrix elements of the relevant operators evaluated below the scale 

~ \J 'Aqcd"^- First we consider an effective theory called SCETi, in which the degrees 
of freedom with p ~ m& are integrated out and the effective theory describes physics below 
the scale m b and above pL . In SCET^ we can construct all the effective operators which are 
gauge invariant under collinear and ultrasoft (usoft) gauge transformations. But a typical 
scale in SCETj is of order m b \ ~ ^/m fc AQ CD , which is still too big for the evaluation of the 
matrix elements of the effective operators. A next step is to go down to the second effective 
theory SCET n , in which all the degrees of freedom with p ~ y mj Aqcd are integrated 
out and a typical scale in SCETn is of order Aq C d and the small parameter becomes of 
order Aq CD /m5. In SCETn, we can evaluate the matrix elements of the operators between 
meson states. The two-step matching looks complicated, but the construction and the power 
counting of operators become manifest. 

Radiative B decays have been previously considered in different theoretical frameworks. 
Ali and Parkhomenko [14] calculated hard spectator contributions to radiative B decays 
in the scheme of the large energy effective theory, and presented their results at next- 
to- leading order and to leading order in AQ CD /m b . In this approach they calculated the 
decay amplitudes in the full theory, and expanded them in powers of 1/E, where E ~ mj. 
Radiative B decays were also considered in the context of the heavy quark mass limit [15, 16], 
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in which the amplitudes in the full theory are expanded in powers of l/m& [17], employing 
the light-cone wave functions of the mesons convoluted with the hard scattering amplitudes. 
The results in Refs. [14-16] show that the radiative B decay amplitudes are factorized at 
next-to-leading order in a s and to leading order in Aqcd/^6- Subleading effects in radiative 
B decays such as the study of the topology of the weak annihilation and the ^-exchanges 
channels [18], and the isospin breaking effects [19] were studied. 

If we apply SCET to radiative B decays, we can extend the proof of the factorization 
theorem to all orders in a s at leading order in Aq CD /m 6 . It is possible because we require 
that the effective operators in SCET n be gauge invariant under collinear and soft gauge 
transformations. The gauge invariance and the reparameterization invariance put a serious 
constraint on the possible forms of the effective operators in SCET [20, 21]. The gauge 
invariance is achieved by introducing the corresponding Wilson lines, and these Wilson lines 
include the effects of collinear and soft gluons to all orders in a s . The Wilson coefficients of 
the effective operators can be calculated to desired accuracy, say, next-to-leading order, or 
next-to-next-to-leading order, but the form of the effective operators remains the same. 

Another advantage in applying SCET is that we can make the power counting explicit 
at each step of the effective theories SCETi and SCETn- In SCETi, we can obtain gauge 
invariant operators in powers of A. When we integrate out an intermediate scale, we need 
not worry about operator mixing in different powers of A since the matching is perturbative. 
As we go down to SCETn, the small parameter now becomes ~ Aqcd/^6 and the power 
counting is explicit. To summarize, our procedure to employ SCET is first to match QCD 
onto SCETi at /i ~ m&, and factorize the usoft-collinear interactions with the field redefin- 
tions. Finally we match SCETi onto SCET n at /i ~ and evaluate the matrix 

elements. 

The organization of the paper is as follows: In Section II, we briefly review the basic 
ingredients of SCET in studying radiative B decays. We explain which types of operators 
are needed in computing the decay amplitudes of radiative B decays, and the outline of the 
procedure in using SCET is sketched. In Section III, all the relevant operators are derived in 
SCETi. There are effective operators without spectator quarks, and if we include the effects 
of spectator quarks, new four-quark operators and four-quark operators with a photon are 
induced. These operators are constructed in a gauge invariant way. In Section IV, the 
operators obtained in SCETi are evolved down to SCET n , and we evaluate all the matrix 
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elements of the operators at leading order in SCET. In Section V, all the contributions are 
combined to produce the decay amplitudes for the radiative B decays. In the final section, 
we summarize all the effects of various operators and a conclusion is presented. 



II. BASICS OF SCET FOR B V-y DECAYS 

In order to describe radiative B decays in SCET, the first step is to construct effective 
operators in SCETi by integrating out the degrees of freedom of order from the full 
theory. The effective operators should be gauge invariant under collinear and usoft gauge 
transformations, and these operators can be systematically expanded in a power series of the 



small parameter A, of order yAQ CD /m b . After we construct all the operators contributing 
to radiative B decays in SCETj, we evolve these operators down to SCETn. In obtaining 



SCETn, we integrate all the off-shell modes of order ^m^KqcTi-, and all the operators in 
SCETn can be expanded in a power seris of the small parameter, which now becomes of 
order Aq CD /m b . At each step of the effective theories, we match theories at the boundary 
requiring that matrix elements of a given operator be the same in both theories. And we go 
down to the scale of interest by using the renormalization group equation. 

We begin with the effective Hamiltonian for B decays in the full QCD, which is given by 

H cS = % £ + C 2 O p 2 + £ QOi), (I) 

V * p=u,c i=3,-,8 

where 

°i = (Pibi)v-A(djPj)v-A, 0\ = (pjb^v-AidiP^v-A, 

3 = (dib^v-Aj^iQjQ^v-A, 04 = (djbi) V -A^2(qiqj)v-A, 

5 = (dibi) V -A^2(qjqj)v+A, O e = {djbi) V -A^{qiqj)v+A, 

Q 1 

Or = -|^vi^(I + 75 )6, 8 = -^d^G^iTaUl + lb )b r (2) 

Here p is an up-type quark u 01 c quark, and d is a down-type quark d (for B — > pj) or s (for 
B — > K*^f) quark. The indices i, j are color indices, F^ v and G^ u are the electomagnetic and 
the chromomagnetic field strength tensors respectively, and Vij are the Cabibbo-Kobayashi- 
Maskawa matrix elements. The sign convention for 07,8 corresponds to negative Cj^ and 
the covariant derivative is defined as = <9 M — igT a A a ^ — ieQA^, where A a ^ (A^) is the 
gluon (photon) field. 



The leading contribution to B — > V'j comes from the operator O7 in the full theory. In 
the effective theory, the leading-contribution comes from the leading-order operator derived 
from O7. In the following computation, we will neglect the radiative corrections of the 
penguin operators since their coefficients are proportional to a s Ci {i — 3, • • • , 6), which are 
numerically small compared to C7. But we include the loop correction from 0\ since a s C\ 
is not negligible compared to Cj. When the penguin operators appear at leading order in 
a s , we include their effects. 

Before we consider the effective operators for radiative B decays in SCET l5 let us intro- 
duce the notations we use in this paper and discuss the power counting of operators. We 
choose the direction of the vector meson as n M , and the direction of the photon as n^. And 
the vector mesons are regarded as massless at leading order in SCET. Therefore the 577(3) 
breaking effects in B — > p7 and B — > K* r y and the SU (2) isospin breaking effects in B — > p°7 
and B~ — > p~7 arise at subleading order. The quark mass effect in SCET was discussed 
in Ref. [22], which can cause leading operators in SCET n . However the matrix elements of 
these operators in radiative B decays at leading order in SCET turn out to vanish. This is 
discussed in Appendix in detail. 

We denote the collinear fields £ and \ as the collinear fields in the n M and the n M directions 
respectively, satisfying 

« = 0,f{ = {,fe = 0,*x = X. (3) 
The collinear gluon field A 1 ^ in the direction can be decomposed as 

A» = -n-A n + A» n± + \n-A n = 0(X°) + 0(A) + 0{\% (4) 

and the power counting of the operators involving collinear fields can be made explicit as 
shown in Ref. [2, 4]. On the other hand, the photon field cannot be decomposed as the 
collinear gluon to facilitate the power counting. The decomposition of the collinear gluon 
field in Eq. (4) is meaningful since quarks and gluons have small fluctuations due to strong 
interaction. If we consider an external photon, the photon field does not have such small 
fluctuations and the momentum of the photon is always of order A , and only the transverse 
part contributes to decay amplitudes. Therefore we do not include the photon field in the 
power counting of operators. Though we assign any power of A to the photon field, it does 
not affect the relative power counting of various operators since all the operators for B — > 
will include a single photon. 



There are two kinds of effective operators in SCETi from different channels, which can 
be classified into the operators without spectator quarks and those with spectator quarks. 
The leading operator for the first type of operators comes from 7 by replacing the quarks 
in the full theory by those in the effective theory and making the resultant operator gauge 
invariant under collinear gauge transformations. We also have the contributions from the 
operators 0\ and 0% by radiative corrections, whose matrix elements are proportional to the 
matrix element of O7. The resultant operators are heavy-light quark bilinear operators with 
a photon field. These operators, after we make them gauge invariant by attaching Wilson 
lines, can produce operators with an additional external gluon. At order g, these operators 
can be obtained by attaching a collinear gluon to O7, or a photon to Og, or a gluon and a 
photon to 0\. 

All these operators without spectator quarks can be obtained by attaching a photon or a 
gluon to fermion lines in the operators 7 , 0\ and 0$. When we attach a photon or a gluon 
to each operator, if any intermediate states have momenta of order A , we integrate out 
these modes to obtain effective operators. If the intermediate states have momenta of order 
\rrib, then these modes cannot be integrated out in SCETi. Instead that channel should be 
considered in SCETn- In order to make the operators gauge invariant under collinear gauge 
transformations, we also have to consider attaching collinear gluons and integrate out all 
the off-shell modes of order A . We compute all the operators in SCETj to order A in order 
to obtain the leading-order result due to the possible enhancement in evolving the operators 
to SCET„. 

The second type of the operators is the operators with spectator quarks. These include 
the four-quark operators of the generic form %Tih ■ £r 2 £, and £Tih ■ xr 2 £. As mentioned in 
Ref. [13], we can use the Fierz transformation to make the operator C,Tih ■ xr 2 £ of the form 
xTih • £r 2 £, and we will consider the four-quark operators only in the form xVih ■ £r 2 £ from 
now on. In these four-quark operators, the spectator quark is a collinear quark x- And we 
can make time-ordered products of these operators with the electromagnetic interaction with 
a collinear quark x and an usoft quark. It is necessary to include the four-quark operators of 
the form x^ih ■ £r 2 £ in SCETi since the time-ordered products of the four-quark operators 
with the electromagnetic interaction contributes to radiative B decays. 

There are also four-quark operators with a photon field of the form q us Y\h ■ £T 2 £ ■ A^, 
where Tj are Dirac matrices. These operators are obtained from the four-quark operators 



6 



in Eq. (1) by attaching a photon to all the quarks except the spectator quark of the B 
meson. In this case, the intermediate state has momentum of m& and can be integrated out 
to produce effective operators. The effects of the operators with a photon attached to a 
spectator quark are already included in the four-quark operators mentioned above. In this 
case the intermediate state has the momentum of order m^A and cannot be integrated out 
in SCETl 

When we evaluate the matrix elements of all the operators in SCETi, we consider the 
matrix elements of the operators without spectator quarks, the four-quark operators with a 
photon, and the time-ordered products of the four-quark operators with the electromagnetic 
interaction of an usoft quark and a collinear quark x- 111 addition, when we consider the 
hard scattering contribution to the form factor for B — > V, we also include the time-ordered 
products of the operators without spectator quarks with the subleading collinear effective 
Lagrangian. And these matrix elements are evolved down to SCET n . By decoupling soft 
gluon interactions, we obtain gauge invariant sets of operators and time-ordered products 
under collinear and soft gauge transformations, and evaluate the matrix elements. 

We can consider other forms of four-quark operators, say, £Tih ■ £^£.4^, in which all the 
light quarks are collinear fields in the n M direction. We can make time-ordered products of 
these operators with the interaction of collinear gluons with an usoft quark and a collinear 
quark £ [21, 23]. If there are no external gluons in the time-ordered products, all the collinear 
gluons should appear as internal gluons, and the final operators correspond to, at higher 
orders in a s , the four-quark operators with a photon of the form q us Tih ■ £r 2 £ • 4. M after we 
perform the time-ordered product with the interaction involving the usoft quark. If there 
are external gluons, these are collinear gluons in the n M direction. In this case, these gluons 
should belong to the final state of a vector meson. Therefore it involves higher Fock space 
for the vector meson state, which we will neglect since the contributions are subleading. 
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III. OPERATORS IN SCETi 



A. Operators without spectator quarks 

In radiative B decays, the leading operator is 7 and the corresponding gauge invariant 
operators in SCETi to first order in A are given by 



(5) 



where A ( f' la,lh " > are the Wilson coefficients, which are operators in SCET. The operator Of 
is the leading operator, and Oj la ' lb ^ are the subleading operators suppressed by A compared 
to Of\ These operators are obtained by attaching a collinear gluon A^ n to the heavy quark 
and integrate out off-shell modes, as shown in Fig. 1. The operators are given as 



Of 

v» 



8^ 



4tt 2 



^^ W )^[WUp n± W}(l + l5 )h, 

e ^^w)[w^ n± w]^(i + l5 )h, 



(6) 



where A 11 is the photon field, and h is the heavy quark field in the heavy quark effective 
theory. The factor W is the Wilson line defined as 

I 



W = ^2 ex P 



-9 



perm 



n-V 



-n ■ A, 



(7) 



where = n ■ Pn M /2 + is the label momentum operator for collinear fields in the n M 
direction. 



On 



FIG. 1: QCD diagram attaching a collinear gluon to the b quark to make a gauge invariant operator. 
A wavy line represents a photon, and a curly line represents a gluon. The momentum of the photon 
(gluon) is outgoing (incoming). 
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When we actually compute the Feynman diagram in Fig. 1, we obtain Of^ and 7 la ^ 
at tree level. However, the operators Oj la ' lb ^ are the most general operators at subleading 
order considering the Dirac structure and the collinear gauge invariance. When we match 
SCET to the full theory at tree level, the Wilson coefficient for 7 lb ^ is zero, but in the 
matching at higher order in a s , a nonzero coefficient A 7 lb ^ can be developed. If we calculate 
decay amplitudes to next-to-leading order accuracy, we need to calculate 4°^ to order a s 
and A 7 la ' lb ^ to order a° s since the operators 7 la ' lb ^ begin with order g and these operators 
will be combined with the effective Lagrangian in the time-ordered products, which are of 
order a s . There are no other subleading operators by including the subleading correction to 
the collinear field £ in Of\ where the fermion field ip in the full theory changes to 

It is because the subleading correction to the collinear field £ is proportional to f, and this 
vanishes when contracted withf£ in the operator 7 °\ 

The Wilson coefficient A 7 °\ to first order in a s , is given as 



4°V) = l-^^ln 2 4 + Jln4-21n4lnn-Pt 



7T 2 i 

+21n 2 n-Pt + 2Li 2 (l-n-Pt) - 21nn-Pt+ +6 , (9) 

where n-V = n- V/rrib, and L^rc) is the dilogarithmic function. The coefficient A^ is the 
same as Cg in Ref. [3]. The remaining Wilson coefficients to zeroth order in a s are given by 

4 la) = i, = o. (io) 

Since the Wilson coefficients are operators, A^ in Eq. (5) means that the Wilson coefficient 
is sandwiched by the operator, that is, it implies that 

4°) Of = ^WfAP(n ■ VW(1 + ls)h = 4 0) (n • p, n ■ q )^Wfifi(l + l5 )h, (11) 

where p (q) is the momentum of the quark £ (the gluon in W). 

We can obtain operators from 0\ and 8 , proportional to 7 by attaching an external 
photon and a gluon in a loop, whose matrix elements are proportional to the matrix element 
of 7 through the radiative corrections. The leading-order contributions from the opera- 
tors 0\ and 8 , which are proportional to Of^ will be included in calculating the decay 
amplitudes in Section IV. 



Now let us consider the operators which consist of a heavy-to-light current with an ex- 
ternal photon and an external gluon, in addition to 0^ a ' lb \ From the operator 0$ in the 
full theory, we can construct such an operator. It is given by 

Os - Q d 4 a) { l a) + QA b) 0?\ (12) 

where Q d is the electric charge of the down-type quark (d or s). Note that there are only 
two independent operators O^' 1 ^ at order A in SCETi. Here the operator without a photon 
from 8 is not included since it does not contribute to radiative B decays, but we comment 
that the leading operator from 8 without a photon field is suppressed by A 2 compared 
to Of\ We can compute the Wilson coefficients A^ a ' lb ^ to any desired accuracy, but at 
next-to-leading order, we need the tree-level Wilson coefficients since the operators begin 
with order g. Matching the operators with the full theory at tree level, as shown in Fig. 2, 
the Wilson coefficient A 8 lb ^ is zero, and QdA^Oj 1 ^ can be written as 

Q d At ] 0? a) = 6 -^W^^n ■ V[W^ n± W](l + 75)/*, (13) 

where Qd is the electric charge of the collinear quark £. The operator in Eq. (13) is obtained 
by attaching a photon to fermion lines in the operator 8 , as shown in Fig. 2. Note that 
the contribution of Fig. 2 (b) is suppressed by A compared to that of Fig. 2 (a), and it is 
neglected here. 

We can also obtain an operator with a photon and a gluon from the operator 0±. We 
cannot make the operator from O2 due to the color structure, and we neglect the contri- 
butions of penguin operators. The full QCD diagram is shown in Fig. 3 and the resultant 




FIG. 2: QCD diagrams attaching a photon to each fermion to the operator Os in order to generate 
the operator with an external photon and an external gluon. The momentum of the photon (gluon) 
is outgoing (incoming). Diagram (a) is of order A compared to o[°\ Diagram (b) is suppressed by 
A 2 and is neglected. 
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operator can be written as [24, 25] 



(14) 



where 



geQu 



8tt 2 



ie^ a p(q a Ai 5 + g£Ai 6 ) 



H — — (2e pi , Q/ 3^g p g£Ai 23 + «e PMQ/ ?<f ^^A^U^l - 7s)- (15) 

Here the momentum of the gluon g is incoming to the loop and the momentum of the photon 
qn is outgoing. The quantities Ai n = Ai n (^ , zi, z 2 ) are functions of the three variables z , 
Z\ and Z2, which are the invariant mass-squared of the fermion-antifermion pair in the loop, 
the gluon momentum squared, and the photon momentum squared, each divided by the 
internal fermion mass m 2 , and they are given by 
(q - Qn) 2 (pi - m b v) 2 



r(l - u), z 1 = — , z 2 



1 2 2\ /' 1 2' z 2 

ml m\ ml m { ml 



0, 



(16) 



where u is the momentum fraction of the d quark in a vector meson. Replacing the quark 
fields with the fields in SCET, it turns out that 1^ starts from terms of order A apart 
from Ai„'s. Therefore we need to expand Ai n to leading order in A, which are given by 
Ai n (zo, 0, 0). In this case, Ai n have the relation 



Ai 23 (z , 0, 0) = Ai 26 (z , 0, 0) = Ai 6 (z , 0, 0) = Ai 5 (z , 0, 0), 



(17) 



and Ai 5 (z , 0, 0) is given by 

4 



Ai 5 (z ,0,0) = 



Li 2 ( 



=) +Li 2 ( 



1 - yjl -A/z + ie J 1 + yJl-4/z + ie' 



=)]+2. (18) 



(a) 




(b) 




FIG. 3: QCD diagrams generating the operator with an external photon and an external gluon 
from Oi, in which a photon and a gluon are attached to the internal fermion line. The gluon 
momentum q is incoming and the photon momentum q^ is outgoing. The same diagrams with the 
penguin operators are neglected. 
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Note that, in Ref. [25], the momenta of the gluon and the photon are both outgoing, while 
the momentum of the gluon is reversed in our case in order to keep a consistent convention 
that a collinear gluon is incoming. Because of this different choice of the directions of the 
momenta, Ai 5 and Ai 2 3 in Eqs. (17) and (18) have opposite signs compared to those in 
Ref. [25]. 

Replacing the fields in Eq. (14) by the fields in the effective theory, and making the 
resultant operator collinear gauge invariant, we obtain the effective operators at order A, 
which are given by 

Or -> Q„4 la) 0? a) + Q u A? b) 0i lb \ (19) 

where Q u is the electric charge of the up-type quark in the fermion loop. To next-to-leading 
order, the Wilson coefficients are given by 

A 0*) = _ A W = lH(j[.Vi,8i), (20) 

where H(n-V\ Sj) = Ai 5 (z — > (1 — n-V' /m^/si, 0, 0) is now an operator with Sj = m 2 /m 2 . 
Compared to h(u,s) in Ref. [16], H(u,s) = —uh(u,s). The operators in Eq. (19) can be 
written as 

Q u A { l a) { l a) + Q u A { l b) 0^ b) = e -^(WH(n.^,s t )^[WUp n± W](l + l5 )h 

-lWH{n ■ V\ Si) [W ] 4> nA _W]^{l + 75 )/i) . (21) 

So far, we have considered the effective operators from the operator Oi, in which a photon 
and a gluon are attached to the internal fermion loop. However, there are other possibilities 
to attach a photon and a gluon to fermions. For example, we can attach a gluon in the 
fermion loop and a photon to external fermions, as shown in Fig. 4. Without the photon, 




FIG. 4: QCD diagrams generating an operator with a gluon and a photon from 0\. The resultant 
operator is suppressed by A 2 compared to , and is neglected. 
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(a) 




(b) 




FIG. 5: QCD diagrams in which a photon attached to the internal fermion line and a gluon is 
attached to one of the external fermions of the operator 0\. These diagrams vanish due to the 
electromagnetic gauge invariance. 

the fermion loop with a gluon attached to it produces an operator [26] 



where rrii is the fermion mass in the loop and q^ is the gluon momentum. For the on-shell 
collinear gluon q 2 = 0, when we replace the quark fields by the fields in the effective theory, 
the operator in Eq. (22) is proportional to £q ■ A n $ ± (l — 7 5 )/i, and it is suppressed by A 3 
compared to Of^ since q ■ A n q^_ ~ 0(A 3 ). When we attach a photon to external fermions 
to make the final effective operators, the power counting does not change, hence they are 
neglected. 

The remaining possibility is to attach a photon to the fermion loop and a gluon is attached 
to external fermions of the operator 0±, as shown in Fig. 5. This operator is zero due to 
the electromagnetic gauge invariance. That is, the fermion loop in Fig. 5 is proportional to 
Eq. (22), replacing the gluon field by the photon field. Since the photon is on-shell, q 2 = 
and q • A — 0, the operator vanishes. 

B. Operators of the form (jTih)(^T 2 

The operators with spectator quarks involve four-quark operators and four-quark oper- 
ators with a photon. The four-quark operators are generically of the form (x^ih)(^2^)- 
The inclusion of these operators in SCETi is necessary because there is an electromagnetic 
interaction of a photon in the n M direction with an usoft quark q us and a collinear quark x- 




(22) 



with 




(23) 
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The time-ordered products of the four-quark operators with the electromagnetic interaction 
correspond to the emission of an energetic photon from the spectator quark. At first sight, 
we can make an operator by attaching a photon to a spectator quark and integrate out the 
intermediate state to make an operator, but the intermediate state has the momentum of 
order m&A, which cannot be integrated out in SCETt. Therefore we construct the four-quark 
operators of the form (xri/i)(£r 2 £) in SCETi and make the time-ordered products of these 
operators with the electromagnetic interaction. Note that, if a photon is attached to exter- 
nal fermions other than the spectator quark, the intermediate state has the momentum of 
order mj,, and it can be integrated out to produce four-quark operators with a photon. This 
will be considered in the next subsection. 

The explicit form of the effective four-quark operators can be obtained from the four- 
quark operators for nonleptonic decays in Ref. [13] with a minor modification. We only have 
to change n M <-> W and x £ m the effective operators in Ref. [13]. These operators are 
given as 



o 1R = 


{(TW) a h a ) v _ A ((?wuw^h) v _ A , 


2 R = 


((TWhK) v _ A {{l d w) a (w^c) P ) V _ A , 


Om = 




4R = 


(tfwhK) y ^ A T,((?w) a (w^%) y _ A 

g 


5R = 


{(x d W) a h a ) v A Y,((?wuw^%) 

g 


Om = 


((x d W) P h a ) v A E{Cewuw^h) v+A 

g 


o lc = 


((^w) P h a ) v _ A ((rwuw^h) v _ A , 


o 2C = 




o 3C = 


E((rw),h a ) v _ A (gtw) a {w^%) y_ A > 

g 


o 4C = 


E {(rw) a h a ) v _ A ((?wuw^%) v _ A . 

g 



(24) 

Here the summation over q goes over to light massless quarks, say, u, d, and s quarks, and 
W is the Wilson line which is given by 

1 

3xp -g- 

L r) 

perm 



w = E ex P 



n ■ Q 



n ■ A w 



(25) 
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b 



^1 



o 



b 



(b) 




(a) 



FIG. 6: QCD diagrams generating the four-quark operators from 0\ and Og. But these operators 
are included in the four-quark operators in Eq. (24). 

where Q 1 = n ■ + Qj_ is the label momentum operator for the collinear fields in the 

n M direction. The operators in Eq. (24) with W and W are collinear gauge invariant. The 
Wilson coefficients of the operators in Eq. (24) are the same as the corresponding operators in 
Ref. [13]. The time-ordered products of these operators with the electromagnetic interaction 
with an usoft quark and a collinear quark x contributes to the radiative B decays. These 
correspond to the contribution from the annihilation and the VF-exchange channels. 

The operators 0\ and Os can induce four-quark operators, as shown in Fig. 6. However, 
these operators are four-quark operators which can be expressed in terms of the four-quark 
operators in Eq. (24). We include these operators in the operators given by Eq. (24), and 
their effects appear in the effective Wilson coefficients of the operators. 

C. Four-quark operators with a photon field 

We can obtain four-quark operators with a photon field by attaching a photon to fermion 
lines in four-quark operators. First these operators can be obtained from the annihilation 
and the IF-exchange channels from the four-quark operators in Eq. (1) by attaching a photon 
on each fermion except the spectator quark. The relevant Feynman diagrams are shown in 




FIG. 7: The four-quark operator with a photon in SCETi arising from the four-quark operators, 
which correspond to annihilation channels or PF-exchange channels. 
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FIG. 8: QCD diagrms to produce the four-quark operators with a photon arising from 0\ and Og. 
These diagrams are subleading. 

Fig. 7. The Feynman diagrams in Fig. 7 give the desired operators in SCETi by integrating 
out the off-shell modes. But these are subleading since the effective operators involve an 
usoft quark. When we go down to SCET n , there are no enhancements and the effect of 
these operators still remains subleading and we neglect these operators at leading order in 
SCET. In Fig. 7, the Feynman diagram in which a photon is attached to a spectator quark 
is missing, but it has been considered in the previous subsection in which we consider the 
four-quark operators of the form (xr 1 /i)(^T 2 ^) and we calculate the time-ordered product 
of this operator with the electromagnetic interaction. This is because the intermediate state 
has momentum of order m b X, and it cannot be integrated out. 

The effective four-quark operators with a photon can also arise from the operators 0\ 
and 0$, as shown in Fig. 8. However, the diagrams in Fig. 8 are subleading due to the 
interaction of a collinear gluon with an ultrasoft quark and a collinear quark. When we 
match this operator to SCETn, there is no enhancement, hence they are still subleading and 
neglected here. 

IV. OPERATORS IN SCETn AND THEIR MATRIX ELEMENTS 
A. Form factors and the leading contributions 

We evaluate the matrix elements of all the relevant operators for radiative B decays, 
which are obtained in SCETi. And we decouple the usoft gluon interactions by redefining 
collinear fields and go down to SCETn. There are various matrix elements we should con- 
sider. First we evaluate matrix elements of the operators without spectator quarks, namely 
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Of^ and the operators induced from 0\ and 8 , of which the matrix elements are pro- 
portional to the matrix element of Of\ Second, we consider the time-ordered products of 
the operators including external gluons without spectator quarks with the interaction of a 
collinear gluon with a collinear quark £ and a usoft quark. These matrix elements contribute 
to the form factor B — > V. As will be seen later, some of the matrix elements are factor- 
ized and they correspond to the contribution from the hard scattering amplitude, which 
is calculable in perturbation theory. The remaining nonfactorizable contributions involve 
endpoint singularities, but these are absorbed into the soft form factor. The time-ordered 
products from other operators, the time-ordered products of four-quark operators with the 
electromagnetic interaction, and the four-quark operators with a photon correspond to hard 
scattering amplitudes with the convolution of the light-cone wave functions of mesons. As 
discussed before, the contributions from the four-quark operators and the four-quark oper- 
ators with a photon are subleading and they will not be included in the final result for the 
decay amplitudes at leading order in SCET. 

In evaluating matrix elements of various contributions, the matrix elements of heavy-to- 
light currents are related to form factors. The form factors can be defined in SCET with 
a fewer number of indepedent form factors due to the symmetry of SCET. In general, the 
form factors for B decays into vector mesons are defined as [27] 



{V(p^)\q^b\B(p B )) 
(V(p,r]*)\qr^b\B(p B )) 



m B + m v 



VuPaPB/3, 



2m v A (q 2 )^q» 
T 

+(m B + m v )Ai(q 2 ) 
rf ■ q 



7/ 



q 2 q 



-Mi 2 > 



Pb+P 



M _ 



m, 



m 



(V(p,OI?^(l + 7fe)6|fl(p B )> = 



m B + m v 
1TM 2 )^rf v p Ba pp 

-tT 2 (q 2 ) \(m 2 B - my)^ - (rf • q)(p B + jf) 



-iT 3 (q 2 )(v* ■ q) 



m, 



u B - m v 

where q = p B — p, my (rj) is the mass (polarization vector) of the vector meson and we use 
the convention e 0123 = —1. Note that the definiton of the form factors are different from 
that defined by Ball and Braun [28] , but if we multiply the right-hand sides of Eq. (26) by 
the factor —i, we obtain the definition in Ref. [28]. 



-(p£+p")l, (26) 
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In SCET, we neglect the masses of the vector mesons, and the polarization vector is 
transverse in this limit. We choose the plane of the transverse polarization defined by 
n ■ r/* = n ■ rj* = 0. In SCET, all the independent form factors defined in Eq. (26) are 
reduced to a few independent nonperturbative functions. Especially, at leading order in 
SCET, there is only one independent form factor (±, and the matrix elements in Eq. (26) 
reduce to [2, 7, 29] 

(V(n,rf)\lW-f&h\~B) = iEC±(E)e^rj* ±v n a n f3 , 
(V(n,ri*)\ZWr^h\B) = 2EC L {E)r t ?, 
{V{n^)\lWa^q u S^h\B) = -Em B C ± (E)e^r ] * ±u n a n fS , 
(Vin^^Wa^q^hlB) = -2iEm B C±(E)vT , (27) 

where the operators are replaced by the operators in SCETn in a gauge invariant form. Note 
that there is another nonperturbative function Q, but since we choose n ■ rj]_ = n ■ rj* = 0, it 
does not appear in Eq. (27) at leading order in SCET. 

From now on, we express all the operators at next-to-leading order accuracy and replace 
the operators of the Wilson coefficients by the Wilson coefficients after we apply the operators 
of the Wilson coefficients inside the relevant operators. The leading-order contribution to 
the radiative B decays comes from the operator Of\ and in SCETn, the gauge invariant 
form under the collinear and soft gauge transformations is given by 

°r 0) = ~^ W W ~^ h - ( 28 ) 
The matrix element of is given by 

(V(n, V± ), 7 (n, e ± )\0^\B) = ^±rn B ( L (E) (2e* ± ■ r)* ± + ie^VWfc) , (29) 

where \B) is an appropriate B meson state depending on the type of the vector meson V . 

There are also radiative corrections at order a s from the operators 0\ and 0§, in which 
there are an external photon and an internal gluon. The matrix elements of these operators 
are proportional to the matrix element of Of\ Denoting the operators induced from 0\ 
and 0$ in as and , their matrix elements are proportional to the matrix element of 
Of\ and can be written as [16, 25] 

(Ofl) = (Of (30) 
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where 



G 1 (s c ) 



104 m t 8 m b 

— In — + 0i (s c ), G 8 = -- In — + 
27 fi 3 /x 



833 20 



162 27 
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— 27T + 



8tt 2 



,3/2 



9 



+ - 48 + 30m - 5tt 2 - 2m 6 - 36C(3) + (36 + 6m - 9tt 2 ) Ins 
9 - 

+ (3 + 6m) ln 2 s + ln 3 s 
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+ 



9s = 



18 + 2tt 2 - 2i?r 3 + (12 - 6tt 2 ) In s + 6m In 2 s + In 3 s 
-9 + 112i7T - 14tt 2 + (182 - 48m) Ins - 126 In 2 s 
+ 



27 

11 2tt 2 2m 
3~' 



* 3 , 



(31) 
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and s c = m 2 /m 2 . There can be subleading operators O^ 10 ' 1 ^ from Oi and 8 in SCET, but 
they start with g 3 . 



B. Contributions to the form factor 

We consider first the effects of the operators o[°' la ' lb \ which are obtained from 7 . We 
calculate the time-ordered products of these operators with the interaction of collinear quarks 
with an usoft quark and a collinear quark £ in SCETV The interaction of collinear and usoft 
gluons starts from order A, but since the propagator of the exchanged gluon is enhanced by 
1/A 2 , we have to include the interaction of collinear and usoft gluons to order A 2 to obtain the 
leading-order result. Among these contributions, the time-ordered products of the operator 
derived from O-j contribute to the form factor, while those from the operators 0\ and 0$ 
contribute to the hard scattering amplitudes. Here we first consider the contributions to the 
form factor. 

The Lagrangian for the interactions of collinear gluons with a collinear quark £ and an 
usoft quark q us is given by [5, 8, 30] 

sy in ■ D n sy in-D n 

4f = l ^^ ^ l Dn / ^ W ^ + h - c -> ( 32 ) 

where 

ig$l = [in ■ D n , ip n ± ], %g$ = \in ■ D n , ip us + % -gn ■ A n \. (33) 
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FIG. 9: Tree-level graphs in SCETi for the spectator contribution to the heavy-to-light form factor. 
The first diagram contributes to Ti^, and the second diagram contributes to 1^2,3,4,5. 



At leading order in SCETi, the time-ordered products, which contribute to the form factor, 
are given as 

T F = J d 4 xT[0?\o)iC^(x)}, Tf = J d 4 xT[O?\0)iC^(x)}, (34) 
Tf = Jd 4 xT[O?\0)ijC^\x)}, T 3 NF = j d*xT[Of\Q)iCf:\x)l 
T, NF = J d 4 xd 4 yT[O?\0)iC^(x)iC^(y)}, T 5 NF = f d 4 xd 4 yT[O?\0)iC^(x)iC^( 



and these are shown in Fig. 9 schematically. After we evaluate the matrix elements of the 
time-ordered products in SCETi, we decouple the collinear-usoft interaction using the field 
redefinitions 

£(0) = yt^ A (o) = Y ^A n Y, Y(x) =Pexp(ig J* dsn ■ A us (ns)), (35) 

in order to go down to SCET n . Matching at /x ~ y/rribA, the usoft fields become soft 
(Y — > S) and the operators are matched onto the operators in SCETn. When we redefine 
the fields using the Wilson line Y, soft gluons decouple in T F . It means that the contributions 
are factorizable to all orders in a s . On the other hand, soft gluons do not decouple in the 
time-ordered products T^ F . Furthermore these contributions have endpoint singularities 
[8]. They are absorbed in the soft form factor (±, as was done in Ref. [8]. 

Among the factorizable time-ordered products, (Tjf ) = (T F ) = 0. The matrix element of 
Tf is given by 

(Tf) = (V n \ J 'd 4 xT[O?\0)i£W(x)]\B) 

= ^^[dn.xfdn.ke™^ 
it N 16tt J J 
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x (V, l\(q s S)(n • ^)^t^(l + 75)^/1 • ^^#(1 + ls)WZ\B) 
= ^^fBhm%(2nl-Z + ie^?Zn°rf) j dr + ^ (36) 

where u (u = 1 — u) is the momentum fraction of the quark (antiquark) inside the meson. 
7£ is the operator which extracts a soft momentum (of order Aqcd) from a soft particle. In 
evaluating Eq. (36), the matrix element involving the B meson can be calculated as 

(0\q s S(n . X )fi(l + >y 5 )S*h\B) 

= J dr+e-^^Ti [^ B (r + )^(l+ 75 )" 

= _lfsms J ^ Tr [i±*W(l + 75)]<^(r + ) 

= if B m B I dr + e- ir ^ 2 ( f>Ur + ), (37) 



where the leading-twist B meson light-cone wave function is defined through the projection 
of the B meson as [27, 31] 

*b{t + ) = - *^[i±*(|t0+( r+ ) +# B (r + )) 75 ]. (38) 

And the momentum-space representation of the vector meson light-cone projection at 
leading order is given by 

M^ = Ml m +Mlp ± , (39) 
where, at leading twist and at leading order in SCET with n • r]\ — n • r]* ± — 



M( = 0, Ml = - l -^Ef ± yi(f) ± (u), (40) 

where 77^ is the transverse polarization of the vector meson. 

Among the nonfactorizable contributions, (T^ F ) = 0. In (T^f), there exist 1/u 2 , and 
l/r + singularities. However, since all endpoint singularities are regulated by Aqcd in the full 
QCD, if all the soft operators are included to cover the endpoint regions, the singularities 
will not arise. We absorb the nonfactorizable contributions to the form factor into the soft 
form factor, given by (±. 
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C. Contributions to hard scattering amplitudes 



There are other time-ordered products of the operators derived from O^s, which con- 
tribute to the nonfactorizable contributions to the radiative decays, which are defined as 

U 1 = J d 4 xT[O?(0)iL^(x)}, U 2 = J d 4 xT[O?(0)ijC^(x)}, (41) 

where the operators O^g are the effective operators derived from Oi >8 at order A. These 
operators are given as 

O? = QJL™0™ + QJi™0™ 

= e -^^WH{n.V,s^\W^p ni _W]{l +1 ,)h 
-IWHin ■ V, Si ) [W^ip n± W]^(l + lb )h) , 

= ^^^^Tpt* ■ HW*ip n ±W](l + 75 )/,, (42) 

which can be obtained from Eqs. (21) and (13). Since the operators start with A, we 
need time-ordered products only with . 
The matrix element of Ui is written as 



{Ui) = (V, 7 | jd*xT[0?(V)iC^{x)}\B) 

a s C F eQ u J dn . x J dn . ke in.m. x /2 



4tt N 167T 

x (V, 7 | (q s S) (n ■ + ■ lWH{n ■ V, Si )^-i^(l + lb )W^\B) 

= C A ^ eQ u f B f±m 2 B (2r]* ± ■ el + it^^e^r^n 13 ) 

OZ7T 1\ v ' 

x [dr + ^±l [duH(m b u,s t )^. (43) 
J r + J u 

If we use the leading-twist light-cone wave function 4>±{u) = 6uu for the vector meson, we 

can replace u by u in the last line of Eq. (43). But in general, if the valence quarks have 

different masses such that the wave function is not symmetric under u <-> u, the present 

form should be kept. 

Note that the matrix element of the time-ordered product with the second operator of 
in Eq. (42), or the operator Oj lb ^ in Eq. (6), vanishes. To prove this, let us consider the 
matrix element of the operator 

( U) = ^ w[w]i p ni _ w] ^ {l + 7s)/l _ (44) 
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The Wilson coefficients, when sandwiched between the heavy quark and the collinear quark 
fields, produce kinematic variables such as n-pi, where p\ is the momentum of the collinear 
quark £. Then the matrix element of the operator Oj lb ^ with the collinear Lagrangian 
is written as 

(V n \ J d*xT[O? b \0)i£U(x)]\B) = / dn-x J dn ■ ke™***<* 

x ((7.5) (n ■ X )_L^(1 + l5 )S^h ■ IW^itil - 7fc)W^|B) (45) 

The matrix element in the last line in Eq. (45) is zero when we use the B mesons and the 
vector meson projections at leading order. Therefore the matrix element of the time-ordered 
products of Oj lb ^ with vanishes. 

Similarly the matrix element of U2 is written as 

(U 2 ) = (V n \ J d'xTp^^iC^^WB) (46) 

= [dn-x [dn-ke^^ 

7T N 16vr J J 

x (V, 7 | (q s S) (n • x)-^^(l + ri&h • {W^M* + 1,)W^\B) 
D. Contributions of the operators (x^ih){^20 

In SCETi, the contributions of the operators (x^ih)(^20 can be calculated through the 
time-ordered product with the electromagnetic current interacting with an usoft quark and 
a collinear quark. The electromagnetic interaction with an usoft quark and a collinear quark 
is given by 

= eQ sp q u jW\, (47) 

where eQ sp is the electric charge of the spectator quark. To be complete, we need to include 
a nonlocal term and the electromagnetic interaction should read 

oM "> zdM^X + T[T g , it™], (48) 
where the operators in the nonlocal term are defined as 

% = x[A n . {Q \ gAw) (0, + 9 M) + «l + 9 M) n . {Q \ gAv / ] \x, 

43 = VusW\$± + 9 M)X + h.c. (49) 
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FIG. 10: Time-ordered products of the four-quark operators with the electromagnetic interaction 
of an usoft quark and a collinear quark \- The circle represents the four-quark operators, and the 
square represents the electromagnetic interaction. 

The operators in Eq. (49) are the same as those denned in Ref. [10] except that n M and n M 
are interchanged. The matching of the operators to the full theory, and the computation 
of the Wilson coefficients have been extensively discussed in Ref. [10]. But the main point 
here is that the electromagnetic interaction at leading order in SCET has the same operator 
structure as given in Eq. (47), but with a different Wilson coefficient including all the effects 
with a nonlocal term. We will show that the contribution of the time-ordered products of 
the four-quark operators with the electromagnetic interaction vanishes at leading order in 
SCET, and we only need the form of the operator in the argument. 

All the contributions from the four-quark operators through the annihilation channels and 
the IF-exchange channels vanish at leading order. The proof is the following. Neglecting 
the color indices, the four-quark operators can be generically written as 

Q i = (xWr 1 h)(gwr 2 w^), (so) 

where r 1; T 2 are Dirac matrices. The time-ordered products with the electromagnetic cur- 
rent, which is schematically shown in Fig. 10, are given by 

M = (V, 7 | J d A xT[Qm^{x)]\B) 

= jdn-xj . qJ^TMB) = 0. (51) 

The matrix element vanishes at leading order in SCET when we take the matrix element of 
the collinear sector £F 2 £ for T 2 — ^(1 ± 75), or 1 ± 75 with the vector meson, they vanish 
using the leading projection for the vector meson. Therefore at leading order in SCET, the 
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annihilation channels and the VF-exchange channels do not contribute to the radiative B 
decays, though there may be subleading contributions. For example, the isopin breaking 
effects in B — > decays are explained by the effects of these operators. 



V. DECAY AMPLITUDES FOR B -> Vj 

Combining all the contributions considered in the previous section, we can write the decay 
amplitudes for radiative B decays. In general, the decay amplitudes can be written as 

Gf em b m 2 B 

: v \^t ± ■ >l± -r it- ■ ■ii ±ll t ±v -itari/3 j 

p=u,c 

x 

/b/i 



r C 7 {m b )Af\^)^C±{E^) + ^fLc i (m b )G i (m b )AM^C±(E^) 

TUB 47T TUB 

/ dudr + T a (fi )J a (u,r + ,fx ,fx)(j)%(r + ,fx)(j) ± (u,fx)], (52) 

where cy = 1 for V = K*, p, and cy = l/y/2 for V = p°, and i — 1,8. Ci(m b ) are the Wilson 
coefficients in the full theory, evaluated at /i — m b , and A^/iq) are the Wilson coefficients 



+ 



in SCETj, evaluated at /i = ym b AQ CD . Gi are the loop corrections, which are given in 
Eq. (31) at order a s . The functions T a (fi ) are the products of the Wilson coefficients in the 
full QCD at fj, — m b , and the Wilson coefficients in SCETi at /i — /io- The functions J a are 
the jet functions which can be obtained in matching SCETi and SCET n . 

The decay amplitudes for radiative B decays at leading order in SCET and at next-to- 
leading order in a s are written as 



AiB^V^ = ^c v J2 ^[(c 7 (m,)4° ) M + ^(C' 1 G 1 + C' 8 G 8 ))(V, 7 |O( 0) |fi) 

V ^ p=u,c " 

+C 7 (m b )(T 1 F ) + C 1 (m 6 )(C/ 1 ) + C 8 (m b )(U 2 )] 

x C 7 (m b )(A 7 '(/j, ) (±{E) + Tra s — m B / dr+ / du _ 

L v m B N m b m B J r + J u 

o JM m b m B J r + J u 1 

+C 8 (m b )(^^( ± (E)G 8 (s p ) 

v 47T 171b 
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| ™ s C F f B f ± ^ f ^ 0+(r + ) /-^>±Hm 



3 A" m b m B J r + 

The Wilson coefficient to next-to-leading order accuracy is given by Eq. (9), and all the 
other Wilson coefficients are 1 at this order because the Wilson coefficients appearing in the 
terms other than the first term in Eq. (53) is proportional to a s . In fact we have to use the 
evolution of the Wilson coefficient A^\/j,) from m b to /i using the renormalization group 
equation because there is a double logarithm. And the evolution of other Wilson coefficients 
should be performed, which is not available yet. However, since the evolution from ji = m b 
to = fio is small, the numerical difference would be small. 

In the final expression of Eq. (53), the first term proportional to Cj is the contribution 
to the form factor B — > V. The nonfactorizable part {T^ F ) has been absorbed in the soft 
form factor (±, and the hard scattering contribution can be calculated. The form factor Fy 
for B — > V as 

(V(n,r] ± )-f(n,e ± )\0 7 \B) = ^m%c v F v 

x (2el • 771 + ie^rj* ± ^n a n p ) , (54) 

where the B — > V form factor Fy is evaluated at q 2 = 0, and in terms of the nonperturbative 
function (±, it is written, at leading order in SCET, as 

F V = A) '{no) (±{E,fi) + -Ka s (n) — m B / dr + / du = , (55) 

m B N m b m B J r + J u 

where nonfactorizable contributions to the form factor from T^ F , which include the singu- 
larities in 1/u 2 and l/r+ are absorbed to the soft form factor (± [8]. This result is consistent 
with the result of Ref. [16], which was obtained at leading order in Aq CD /m fe in the heavy 
quark limit. Therefore we conclude that the leading-order result in SCET corresponds to 
the leading-order result in the heavy quark mass limit, as it should be. 

In Ref. [16], they consider the scale dependence of the matrix element of Oj by including 
the running of the product of the b quark mass and form factor as 

(m b F v ) M = (m b F v ) [m b ] (l + ^8C F In ^) , (56) 

v 47T /! ' 

but in our formalism, we use the running b quark mass at /i = m b and the Wilson coefficient 
A^ evaluated at /x , and the scaling of Fy is manifestly expressed in Eq. (55). 
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VI. CONCLUSION 



We can apply SCET to the radiative B decays such as B — > K*^y or B — > p7, and 
the organization of the relevant operators can be systematically achieved using the power 
counting in SCET. Employing the two-step matching from SCETi to SCETn, off-shell modes 
are integrated out successively and we can obtain gauge invariant operators under collinear 
and soft gauge transformations in powers of A. The scaling of A changes from 
in SCETi to Aqcd/^6 in SCETn, and the power counting of the operators can be achieved 
including the time-ordered products of the operators with the Lagrangian. The matrix 
elements of the operators can be evaluated in SCET n including the effects of the evolution 
of the operators. We considered four-quark operators through the annihilation channels and 
the W exchange channels in SCET, and all these effects turn out to be at subleading order 
by our power counting method. 

All the hard scattering amplitudes are shown to be factorized and the decay amplitudes 
can be written as the convolution of the hard scattering amplitudes with the light-cone wave 
functions of the mesons. For the contribution to the form factor, we can categorize the time- 
ordered products of Oj°' la ^ with the interaction of a collinear gluon with a collinear quark 
£ and a soft quark into factorizable and nonfactorizable contributions. The factorizable 
contributions are calculable as a convolution of the hard scattering amplitude and light- 
cone wave functions of mesons. In the nonfactorizable contributions to the form factor, 
the effects of soft gluons are not decoupled and they are purely nonperturbative effects, and 
they have endpoint l/u 2 , l/rj. singularities. These nonfactorizable contributions to the form 
factor are absorbed into the definition of the nonperturbative form factor 

We have explicitly shown that the decay amplitudes are factorized to all orders in a s at 
leading order in SCET. This is an extension of the factorization theorem proved at next- 
to-leading order in the heavy quark limit [16]. It is also shown that the decay amplitudes 
at leading order in SCET coincide with the amplitudes obtained in the heavy quark limit 
at leading order in AQCD/^b- This should be true because the two limits at leading order 
are the same, and this correspondence was also shown in nonleptonic decays into two light 
mesons [12, 13]. We have calculated the decay amplitudes at next-to-leading order. But 
this does not mean that the factorization theorem is valid only to this order. The operators 
obtained in SCET are gauge invariant by attaching the Wilson lines, and these include 
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the effects of gluons to all orders. If we analyze radiative B decays at next-to-next-to- 
leading order, the only changes occur in the Wilson coefficients and we have to include more 
time-ordered products, the form of the operators does not change. And the factorization 
properties remain intact. Therefore the factorization theorem works to all orders in a s at 
leading order in SCET, and we have evaluated the Wilson coefficients at next-to-leading 
order in this paper. 

The contributions of the annihilation channels and the l^-exchange channels are sub- 
leading. This point was also pointed out in the heavy quark mass limit [16]. Therefore the 
isospin breaking effects due to the different contents of the spectator quarks in the B me- 
son appear only at subleading order in AQcr>/ m b i n the heavy quark limit [19]. The S77(3) 
breaking effects which arise in the difference between the decay rates B — > 7T*7 and B — > p7 
also appear at subleading order. It would be challenging to consider radiative B decays to 
order Aq CD in SCET. In order to probe subleading effects such as the isospin breaking ef- 
fects and the 577(3) flavor breaking effects, we have to devise a systematic method to derive 
gauge-invariant subleading operators and the time-ordered products at a given order. Re- 
cent discussion [21, 22] on the field redefinitions to make the gauge invariance explicit to all 
orders, and the effects of the quark masses would help construct the operators. 
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APPENDIX: QUARK MASS EFFECTS IN RADIATIVE B DECAYS 



In Ref. [22], quark masses are included in the SCET Lagrangian. If the strange quark 
mass scales as Aqcd, the strange quark mass terms are suppressed in SCETi, but are leading 
order in SCETn- The leading-order Lagrangian is given as 

C = l[in -D+(f ± + g^ nA )wJ—WUf ± + ^n±)]|e, (A-l) 
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FIG. 11: Feynman diagram for the time-ordered product T m from C rn . The square represents 
and the cross represents C m . The momentum q 2 of the gluon is incoming to C m . 

and the mass terms are given by 

Cm = mZ[(f± + 9^n±) , W^-wf-i - m^wJ—W^h. (A.2) 
lv / n -V J 2 n -V 2 

We will consider only the case of the strange quark since the up and the down quark masses 
are very small. As discussed in Ref. [22], if the strange quark mass scales as Aqcdj C m is 
suppressed. However, when we go down to SCETn, there is an enhancement and we have 
to consider the operators to order A in SCETj to obtain a leading-order result for radiative 
B decays. Therefore we include the effect of the first term in C m . 

The additional time-ordered product due to the mass terms at leading order in SCETi is 
given by 

T m = J d A xd 4 yT[O?\0),tC m (x),tC^(y)}. (A.3) 

The Feynman diagram for T m is given by Fig. 10. The time-ordered product in which there is 
no gluon from the vertex given by C m becomes zero when we match the operator to SCETn 
because it is proportional to pj_. The time-ordered product is written as 

Tm = ^±g 2 m^ [ d4xd4<l1 d *y d ** 2 c -i<n-*-in-(x-v)j n -x c -ik-«j n * 



4tt 2 " N J (2tt) 4 (2tt) 4 

xliil^J— - _ , l — , 7iJ - l5)hj\q US!j ^i, (A.4) 
^ n-qi n- {q 1 + q 2 ) > q{ qi ,J 



where the momenta are specified in Fig. 11. Here i, j are color indices, and we omit the 
Wilson lines, but after the calculation, we can always make the operator gauge-invariant. 

In going down to SCET n , after we decouple soft gluons as was done in Eq. (35), we 
can evaluate the matrix element of the resultant operator. However, the basic form of the 
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operator to be evaluated is given by 

?i7±^(l- 7 5 )V^7^ = ^7m(1 + 7 5 )A'-^(1-7 5 )^ 

+^,7m(i + 7 5 )A- • g SJ Y(i + 75)6 
= -1.(1-75)^.^.(1 + 75^. 

+ ^7m(1 + 75)6 • ^"(l + 7s)A', ( A -5) 

where the Wilson lines and the momentum factors are omitted for simplicity. When we 
take the matrix elements of the final operators in Eq. (A. 5), they vanish at leading order in 
SCET. The quark mass effect may appear at higher orders in a s or at subleading order, but 
at the order we consider, there is no quark mass effect. 
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